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Miclielle Mills Strout, Larry Carter, Jeanne Ferrante, Barbara Kreaseck 
February 2004 International Journal of High Performance Computing Applications, 

Volume 18 Issue 1 
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In modern computers, a program's data locality can affect performance significantly. This 
paper details full sparse tiling, a run-time reordering transformation tliat improves the data 
locality for stationary iterative methods such as Gauss-Seidel operating on sparse matrices. 
In scientific applications such as finite element analysis, these iterative methods dominate 
the execution time. Full sparse tiling chooses a permutation of the rows and columns of the 
sparse matrix, and then an o ... 

Keywords: computer architecture, data locality, irregular grids, iterative alogorithms, 
sparse matrix, static and dynamic analysis, tiling 
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We consider the formulation of the Schwarz alternating method for a new class of elliptic 
cubic spline collocation discretization schemes. The convergence of the method is studied 
using Jacobi and Gauss-Seidel iterative methods for implementing the interaction among 
subdomains. The Schwarz Cubic Spline Collocation (SCSC) method is formulated for 
hypercube architectures and implemented on the NCUBE (128 processors) machine. The 
performance and convergence of the hypercube SCSC algorithm is s ... 

1 1 Reliable computations and their applications (RCA): On considering an interval 
constraint solvin g algorithm as a free-steering nonlinear Gauss-Seidel procedure 
Frederic Goualard 

March 2005 Proceedings of the 2005 ACM symposium on Applied computing 
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We show that a classical interval constraint propagation algorithm enforcing box 
consistency may be interpreted as a free-steering nonlinear Gauss-Seidel procedure. This 
suggests that the choice of a transversal in the incidence matrix associated witli the 
problem to solve is paramount to the efficiency of the algorithm. We present experimental 
evidences that it is indeed so, and we suggest an heuristics to compute good transversals. 
The improved interval constraint algorithm is compared with a c ... 

Keywords: branch-and-prune method, constraint, nonlinear system 
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Inverse modeling is an application for adaptive filters that has found extensive use in many 
engineering disciplines. In this paper, we consider the problem of finding reverse models in 
the area of channel equalization, and adaptive control systems. First, the problem is 
formulated in a general setting as a standard least squares problem. With this, the inverse 
model can be found using any one of the many well established least squares methods. One 
such method is the classical Gauss-Seidel metho ... 

Keywords: Gauss-Seidel method, SOR method, adaptive control system, channel 
equalizer, linear acceleration 



1 3 Parallel trian g ular decompositions of an oil refinin g simulation Q 
Xiaodong Zhang 

August 1993 Proceedings of the 7th international conference on Supercomputing 
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One important process In oil refining is to separate the crude oil into various oil products. 
This process Is called distillation. In designing a complex distillation column, a large 
computer simulation is conducted. This paper presents our experience with parallelizing an 
oil refining simulation application that computes the composition of the various oil products 
in designed refining columns operated under a given set of conditions. Mathematical models 
for the simulation ... 

14 Virtual path bandwidth allocation in multiuser networks Q 
Aurel A. Lazar, Ariel Orda^ Dimitrios E. Pendarakis 

December 1997 IEEE/ ACM Transactions on Networking (TON), volume 5 issue 6 
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Keywords: Nash equilibrium, bandwidth allocation, distributed algorithms, game theory, 
network control, virtual path 
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Claudia Leopold 

March 2002 Proceedings of the 2002 ACM symposium on Applied computing 
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Full text available:^ 



Iterative solvers such as the Jacobi and Gauss-Seidel relaxation methods are important, but 
time-consuming building blocks of many scientific and engineering applications. The 
performance problems are largely due to cache misses, and can be reduced by tiling the 
codes. Whereas previous research has shown the usefulness of tiling by experimentally 
comparing the run times of tiled and original codes, it did not tackle the question as to 
whether further improvements are possible. In this paper, we ... 

Keywords: data locality, lower bounds, relaxation methods, tiling 
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Real-world objects are usually composed of a number of different materials that often show 
subtle changes even within a single material. Photorealistic rendering of such objects 
requires accurate measurements of the reflection properties of each material, as well as the 
spatially varying effects. We present an image-based measuring method that robustly 
detects the different materials of real objects and fits an average bidirectional reflectance 
distribution function (BRDF) to each of them. In or ... 

Keywords: BRDF measurement, normal map acquisition, photometric stereo, shape from 
shading, spatially varying BRDFs 
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The hunnan hand is a masterpiece of mechanical complexity, able to perform fine motor 
manipulations and powerful work alike. Designing an animatable human hand model that 
features the abilities of the archetype created by Nature requires a great deal of anatomical 
detail to be modeled. In this paper, we present a human hand model with underlying 
anatomical structure. Animation of the hand model is controlled by muscle contraction 
values. We employ a physically based hybrid muscle model to convert ... 

18 Interactive global illunnination: Interactive global illumination using selective photon Q 
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Kirill Dmitriev, Stefan Brabec, Karol MyszkowskI, Hans-Peter Seidel 

July 2002 Proceedings of the 13th Eurographics woricshop on Rendering EGRW '02 
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We present a method for interactive global illumination computation which is embedded in 
the framework of Quasl-Monte Carlo photon tracing and density estimation techniques. The 
method exploits temporal coherence of illumination by tracing photons selectively to the 
scene regions that require illumination update. Such regions are Identified with a high 
probability by a small number of the pilot photons. Based on the pilot photons which require 
updating, the remaining photons with similar paths I ... 

1 9 Shadin g and shaders: Efficient renderin g of spatial bi-directional reflectance | 
distribution functions 

David K. McAllister, Anselmo Lastra, Wolfgang Heidrich 

September 2002 Proceedings of the ACM SIGGRAPH/EUROGRAPHICS conference on 
Graphics hardware 
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We propose texture maps that contain at each texel all the parameters of a Lafortune 
representation BRDF as a compact, but quite general surface appearance representation. 
We describe a method for rendering such surfaces rapidly on current graphics hardware and 
demonstrate the method with real, measured surfaces and hand-painted surfaces, We also 
propose a method of rendering such spatial bi-directional reflectance distribution functions 
using prefiltered environment maps. Only one set of maps is ... 
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Keywords: graphics hardware, reflectance & shading models, rendering hardware, texture 
mapping 
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This paper presents several numerical methods which may be used to obtain the stationary 
probability vectors of Markovian models. An example of a nearly decomposable system is 
considered, and the results obtained by the different methods examined. A post mortem 
reveals why standard techniques often fall to yield the correct results. Finally, a means of 
estimating the error inherent in the decomposition of certain models is presented. 

Keywords: Markov models, near-decomposability, numerical techniques, simultaneous 
iteration 
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